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Short-range correlations in Ag-Pd alloys are investigated by analyzing the ab initio total energy 
of fee based random AgcPdi_c. Since the information on the atomic interactions is incorporated in 
the energetics of alloys it is possible with a suitable model, Bethe-Peierls-Weiss model is used in the 
present work, to invert the problem, i.e. to obtain information on the short-range correlation from 
the total energy of a random system. As an example we demonstrate how site correlations can be 
extracted from random alloy data. Bethe-Peierls-Weiss model predicts positive first neighbor corre- 
lator and mixing energy for substitutional face centered cubic (fee) Ag-Pd alloys at low temperature 
which can be related to the optimal structures of Ago.sPdo.s. 
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I. INTRODUCTION 

The low-temperature short-range order of Ag-Pd 
has attracted theoretical research for several decades. 
Results supporting phase separation^-^ as well as 
orderin^.^-*^ have been reported. Due to the low transition 
temperature predicted for the disorder-order transition 
there are no direct experimental observation concerning 
this matter. 

In the present work, we reanalyze the Ag-Pd binary al- 
loy system using a spin lattice model beyond the simplest 
mean field (Weiss) model. We employ Bethe-Peierls- 
Weiss (BPW) modcl^S, which, in contrast to Weiss 
model, incorporates non-vanishing correlators of neigh- 
bor sites. Therefore it is suitable for extracting ordering 
information on the binary system in hand. 



II. BPW MODEL 

The lattice structure of binary alloys can be modeled 
by spin (Ising) lattices identifying spin states ±1 of a 
particular site with occupation of a given atom, say A 
or B. For spin model calculations one has to specify the 
interaction energies between the sites. This can be done, 
for instance, by fitting the calculated lattice quantities to 
the corresponding observed or simulated ones. 

In the following we model Ag-Pd by an Ising-lattice 
with N sites each associated with spin — ±1, i — 
1,...N. We identify the spin state cr^ = +1 with the 
occupation of the site i by a Ag atom and (Tj = — 1 with 
the occupation of the site z by a Pd atom. The average 
concentration of the whole lattice is fixed by requiring 
that there is sites with ai = +1 and 7V_ sites with 
ai ~ —1 (iV+ + iV_ = N). The average concentration of 
Ag atoms is then c = N+/N and the expectation value 
over all lattice sites is 



sum IS 



N 



N 



= 2c- 1. 
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where the sum is over all spin configurations with con- 
straint J2i '^i = (2c— 1)A^. The total energy [E) is writ- 
ten as the sum of pairwise interaction energies 



(3) 



An important quantity measuring the ordering of spins 
is the correlator 



[aiai 



(4) 



Here gij G [—1,1], gi,j < (> 0) corresponds to spins 
at sites i and j tending to align antiparallel (parallel), 
whereas the case gij = corresponds to completely ran- 
dom alignment. 

In practise to calculate the state sum in Eq. 1^ one 
has to make further approximation. In the following we 
use the nearest neighbor (NN) approximation, where only 
the closest atom sites are included in the energy sum ^ . 
While Ising models include only two-site interactions, our 
task requires to include the multi-site interactions to the 
model. Formally this is done by allowing the Ising model 
parameters to depend on the concentration c, so that 
many-site interactions are effectively and on the average 
taken into account. 

Thus in NN-models only the nearest neighbor sites con- 
tribute to the energy and for a fixed concentration c only 
effective energy for a site is given. For that purpose we 
define the pair interaction energy as 



e(c) + Ae(c) aiaj + ^eic) {a^ + aj) 




Considering only two-particle interactions the state 



i,j are NN 
otherwise 
(5) 

Suppose now, that the site i of the Ising-lattice has ui 
nearest neighbors j G A/i. This number Ui is the coor- 
dination number of the site i. Then the energy of the 
whole system is 



i=l jeAfi 
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where 



E, 



(7) 



is the effective energy of the site i. In this sum there is 
Vi addends and reads for ^ as 



E, 



1 



1 



(8) 



To go beyond mean field models we treat the system by 
BPW model^*^. In BPW model each site i interacts with 
its v nearest neighboraii il, . . . iiy, which in part interact 
with their NN's other that i itself, i.e. with v — 1 sites. 
The interaction energy of neighbors of i are included ex- 
actly, whereas the interaction of the neighbors with their 
neighbors is calculated using mean field. Thus, instead 
of Eq. d?]) the energy of a site i is 



E, 



E 

i^e + -ei/CTj + (AecTj + -e)[aii H 

+fc(i^- l)[i/e+ (Ae(cr) + ^e){aa + 



a,,) (9) 



where fc is a combinatorial parameter used to obtain spe- 
cific thermodynamic quantities by differentiation and the 
choice k = 1 corresponds to the normal BPW model. For 
Ag-Pd the pair interactions just at the first two coordi- 
nation shells are needed to qualitatively understand the 
ordering energy^. 

For v NNs of site i we denote the number of aij — 
+1 by 71+ and the number of atj = —1 by n_, so that 
+ rt_ = ly and 

Ei = E{at,n+) = ly e+ ^eiyai + {Aeai + ^€){n+ - n^) 



+ kiiy-l) 



-iy[-e- Ae]ai 



AecTi -I- -e + k{v- 1) <^ Ae{a) + -I 



1 



2n+.(10) 



When <C -/V+ and v <^ N we may approximate that 
configuration of the v neighbors of a fixed site i is effec- 
tively independent on all other sites. This is the main 
idea of Bethe-Peierls approximatior^i^. However, BPW- 
approximation is widely used to model short-range cor- 
relation in various statistical system a^'^'^"^ . The validity 
of this approximation is however not clear, as there is 
no known error estimation method; PBW model has an 



uncontrolled error. Anyway, with this assumption of ef- 
fective independence of neighbors of separate sites leads 
to partition function 



N 



N 



(11) 

where /3 = {k^T) ^, k^ is the Boltzmann constant and 
T temperature. 

The effective one site (neighbor) partition function for 
spin a can be calculated using grand canonical ensemble 
of its neighbors. We write 



Zii'y) = E 



The chemical potential of a site /j,' has to be related to 
the overall concentration c, as is still fixed by overall 
concentration condition c — N^/N. 

Here the partition function Z is a grand canonical par- 
tition function with respect to the number of cr = +1 sites 
of the system but canonical one with respect to the total 
number of sites. That is, the partition function is related 
to grand potential 



n 



In Z. 



Thus when we later turn to use Helmholtz free energy 
(i.e fixed concentration) we have to make the appropriate 
Legendre transformation. 

The thermodynamical quantities are to be calculated 
from the logarithm of the partition function, i.e. 



InZ 



In 



N 



1 



Nl3[Eo + cAE{+l) + (1 - c)AS(-l)] 



+ |i/c In 

+ (1 - c) In 



1 + ze 



1 + ze' 



(13) 



with 



En = V 



+ k{v-l){e-Ae{a)--l{l-{a)) 



i^^e — iy*-e + k{v — 1) I -e — Ae 



AE{a) 
E{a) 



AE a 



^2 



Ae]<7, 



Ae cr + + k{v - 1) | Ae(cr) + 
Aecr + vj^e + k{v - l)Ae{a) 
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where z = e'' is the one-site fugacity 



III. RESULTS 



(o-)2 + (1 - ((7)2)e-'3(-E+-^-) - {a) 



1 + M 



(14) 



E± — E{±1), and i^* = 1 + fc(i^ — 1). The internal energy 
reads 



U 



cvE+ {l-c)vE. 



N 



(15) 



Now we turn to use fixed concentration, whence the 
Helmholtz free energy is given by 

which can be expressed as energy density. The formula 
for the entropy can be given after that in a standard way: 

TS = U - F = U +-AtlZ --N+u\nz 

P P 

Nv 

= T5o + — [c/+ + (l-c)/_ -clnz], (17) 



where 



pEi 



+ ln[l + ze-''^±]. (18) 



1 + z-ie^^± 
Further, the entropy per site is 

s ~ So + ks V [c/+ + (1 — c)/_] — ksv c\nz (19) 

and the Helmholtz free energy per site reads as 

/ = -^lnZ+i.clnz. 

The mixing energy per atom can be constructed straight- 
forwardly as 

Umix{c) = u{c) ~ CU{1) - (1 - c) w(0). 

By differentiation of (jl3p with respect to energy pa- 
rameters e and Ae we find 



Next step is to determine the BPW parameters e, e 
and Ae by fitting the internal energy to the calculated 
total energy of AgcPdi_c alloys. The reference data we 
are using are the total energies from Ref. [l^ calculated for 
AgcPdi_c using the exact muffin-tin orbitals metho d^^i^^ 
within the coherent potential approximation^^'—, i.e. cor- 
responding to the mean field approximation (Table |T]) . 
Because the ab initio energies in Table U refer to the non- 
correlated Ag-Pd system at K temperature we have 
to use in the fitting procedure the non-correlated BPW 
energy obtained by differentiation with respect to k or 
equivalently from Eq. © as 

N 
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assuming that {<7i(jj) = {ai){aj) = (cr)^ (i ^ j). We 
find the non-correlated internal energy Unc to be up to 
normalization same as in the Weiss model, 



(21) 



This can be expressed as a function of concentration as 



U„c{c) 



1 



(22) 



where the parameters are related by 



Ae ^ ^A 



(Aec + ec), 
1 1 



(23) 



Ae 



"nc(c)[Ry] 



0.00 
0.05 
0.10 
0.20 
0.40 
0.50 
0.60 
0.70 
0.80 
0.90 
1.00 



10084.720387 
10111.731714 
■10138.743066 
■10192.765913 
■10300.812124 
■10354.835440 
■10408.858240 
■10462.880386 
■10516.901915 
■10570.922942 
■10624.943741 



9i,ii = [1 - (ct) 



1 



1 



l + z-ie/3S+ 1 + z-^eP^- 



(20) 



Thus BPW model clearly allows correlation between a 
site and its neighbors. 



TABLE I: The ab initio internal energy of AgcPdi_c at K 
as obtained in Ref. Iisl. 



To fit the Eq. ((22|) to the numerical data shown in Ta- 
ble |I] is not an unambiguous procedure. Several justified 
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strategies can be introduced to define the concentration 
dependence of the parameters Cc, Cc, and Aec- We study 
four different fitting schemes and compare their results 
with the data obtained from impurity calculations for Ag- 
Pd system. In three first fitting schemes we take Cc and 
Cc to be constants whereas different boundary conditions 
for Aec are introduced. 

i) The choice Aec(l) — Aec(O) = leads to the cor- 
rect linear form of the total energy in the limit of 
noninteracting atoms. This leads to tenth degree 
interpolation polynomial for Aec(c). 

ii) Since the term c^Aedc) vanish at c = the above 
constraint can be lifted to Aec(l) = 0. This leads to 
ninth degree interpolation polynomial for Aec(c). 

iii) Fitting of eight degree polynomial to the data in 
Table |l] and interpreting the two terms of lowest 
degree as ec and ecC, now neither AedO) = nor 
Aee(l) - 0. 

iv) Redlich-Kister (RK) parameterization^^ for Eq. 



e[Ry] Ae[/iRy] 



Unc{c) = -VI^4<^RK + £rkc+ c{l - c)AeiiK{c)]. (24) 



Part of the linear term in ([2^ is included to the 
RK-interaction term. Note, that here e^x and Irk 
are constants. Thus Cc = ^rk, Ec — ^rk + Asrk 
and Aec(c) = —Acrk- The RK-parameterization 
is symmetric in the exchange of PdoAg. 

To decide which one of the above fitting procedures de- 
scribes best the energy of Ag-Pd alloy we compare the 
parameters e(c), e(c), and Ae(c) at both ends (c = 0, 
c= 1) with the corresponding calculated data of Ag-Pd. 
The reference parameters are shown in Table |IT1 





c = 


c = 1 


Ecomp [Ry] 


-143.8171426 


-143.8171231 


Ecomp [Ry] 


-3.7515816 


-3.7516004 


Aecomp l/J-'Ry] 


0.0 


60.417 



TABLE II: Energy parameters e, e and Ae obtained from 
Refs. HHandlll. 

It turns out that the fitting procedure (iii) gives the 
best overall agreement with the ab initio impurity data. 
The BPW parameters obtained from the fit (iii) are 
shown in Table IIIII Using the BPW parameters various 
thermodynamic and statistical quantities can be calcu- 
lated for AgcPdi_c. Figures [T] and [5] show the mixing 
energy and the neighbor correlator at different tempera- 
tures. 

At low temperatures the mixing energy and the neigh- 
bor correlator are positive suggesting that the BPW mul- 
tisite interactions drive substitutional fee Ag-Pd to favor 
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817139 
817143 
817147 
817154 
817159 
817159 
817153 
817144 
817135 
817128 
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■3.75158336 
■3.75159123 
■3.75159910 
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■3.75162327 
■3.75162367 
■3.75161104 
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■3.75155107 



-5.555558 

-9.490742 

-13.425927 

-21.296297 

-25.513510 

-25.712963 

-19.396219 

-10.747355 

-2.199074 

5.178326 

10.587963 



TABLE III: Energy parameters e, e, and Ae as a function of 
concentration according to the fitting scheme (iii). 



an atom and its neighbors to be of the same type. This 
suggests phase separation or segregation tendency for Ag- 
Pd alloys within substitutional fee structures at low tem- 
peratures in agreement with the results of Refs. [l] andH. 
However, as the recent investigations of Delczeg-Cirjak 
et al.'^'^ show the ground state structure of Ago.sPdo.s is 
not a substitutional fee type structure but the Lli type 
structure with c/a larger than its ideal value. Here c and 
a are the lattice parameters in hexagonal representation. 
In conventional cubic coordinate system c axis is along 
the [111] direction and the lattice parameter a is in the 
(111) plane. For ideal fee structure c/a = ^/6, where c is 
the diagonal of the conventional cube and a is the nearest 
neighbor distance. Since the Lli structure is composed 
of alternating Ag and Pd (111) layers the increasing of 
c/a compared to the ideal value means, that the near- 
est neighbor distance of unlike pairs becomes larger than 
that of like pairs. This shifts more weight on the like-pair 
interaction than on the mixed-pair interaction. There- 
fore, the results of BPW model are consistent with the 
predicted tendency of Ago.sPdo.s at lower temperatures 
to order in Lli structure with increased c/a. 

At medium temperatures (50 - 100 K) the mixing en- 
ergy is negative but the maximum at about c = 0.5 sug- 
gests a slight tendency of phase separation. At high tem- 
peratures the BPW mixing energy approaches the ex- 
perimental mixing enthalpy^! which has its minimum at 
about c = 0.6. 



IV. SUMMARY 

We have used the Bethe-Peierls- Weiss model to inves- 
tigate the effect of multisite interactions on the ordering 
in Ag-Pd alloys. The mean field ab initio data has been 
used to determine the parameters of the BPW model. 
The BPW mixing energy and neighbor correlator for a 
substitutional fee structure of Ag-Pd are positive at low 
temperatures supporting the stability of the Lli struc- 
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ture of Ago sPdo 5 with elongation along the hexagonal FIG. 1: (Color online) Helmholtz free mixing energy per site 
[111] axis. frm. of Ag,Pdi_, at T = IK, 20K, 50K, lOOK, 200K (from 

up to down) . 
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FIG. 2: (Color online) Neighbor correlators of AgcPdi_c at 
various temperatures T = IK, 20K, 50K, lOOK, 200K (from 
up to down). 



